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Abstract —We provide convex necessary and sufficient conditions for 
the robust stability of linear positively dominated systems. In particular 
we show that the structured singular value is always equal to its convex 
upper bound for nonnegative matrices and we use this result to derive 
necessary and sufficient Linear Matrix Inequality (LMI) conditions for 
robust stability that involve only the system’s static gain. We show how 
this approach can be applied to test the robust stability of the Foschini- 
Miljanic algorithm for power control in wireless networks in presence of 
uncertain interference. 


I. Introduction 

A dynamical system is said to be positive if, for every nonnegative 
input, the output trajectory remains nonnegative for all time m 
Positive systems have received increasing attention in recent years, 
not only for their practical relevance but also for their interesting theo¬ 
retical features. The study of positive dynamical systems is intimately 
related to the spectral theory of positive and nonnegative matrices 
that dates hack to 1907 with Perron Q and 1912 with Frobenius (3. 
In fact positive linear systems present a simple dynamical behavior 
that evolves along the Perron-Frobenus eigenvector of the matrix that 
describes their dynamics. More recently positive systems have been of 
interest in the contest of biology where several complex phenomena 
can be modeled naturally using compartimental systems @1, whose 
dynamics is inherently positive. The simple dynamical behavior of 
positive systems simplifies greatly their analysis leading to strong 
and scalable results. For example, for these systems. The design of 
static output feedback controllers can be done with linear program¬ 
ming m, the Kalman-Yakubovic-Popov (KYP) Lemma simplifies 
and a diagonal Lyapunov matrix is enough prove contractiveness m- 
This results leads to the design of optimal Hoo distributed controllers 
using convex optimization. In m the authors show that the KYP 
lemma is equivalent to a linear program and therefore both Hoo 
analysis and synthesis can be extended to large scale systems. In H 
the authors introduce a class of systems called positively dominated 
which is a proper superset of the class of positive systems and 
yet maintains most of their attractive properties despite showing a 
richer dynamical behavior. Second order oscillators, for example, are 
positively dominated provided they have sufficiently large damping. 
In this paper we mostly focus on the robustness analysis for positively 
dominated systems. The problem of robustness has been tackled for 
positive systems by El and Go) in the Li and Loo gain setting 
where the authors, exploiting linear copositive Lyapunov functions, 
provide tractable necessary and sufficient conditions for robustness 
analysis and controller synthesis, and in GD where the authors 
consider norm-bounded static unstructured perturbations. In this work 
we focus on the structured singular value setting where we consider 
structured mixed real LTI uncertainties. The structured singular value, 
/i, was introduced by Doyle in GD and provides necessary and 
sufficient conditions for robust stability of LTI systems in presence of 
structured norm-bounded LTI uncertainty. However such conditions 
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are, in general, NP hard to verify. In this work we show that for 
positive and positively dominated linear systems these conditions 
are convex and thus easily computable. Early work showing that 
the structured singular value is tractable for positive systems in the 
specific case of diagonal perturbations can be found in G3l- The main 
contributions of the paper can be summarized as following: 

• We show that, for positively dominated systems, the worst-case 
structured LTI perturbations are alway static real and nonnega¬ 
tive. Furthermore we show that repeated scalar-times-identity 
uncertainty blocks can be treated as full block uncertainties 
without loss of generality. 

• We show that for nonnegative matrices the structured singular 
value is equal to its convex upper bound for general block 
structures involving mixed real full and repeated scalar blocks. 

• We show that for positively dominated systems robust stability 
with respect to a structured mixed real perturbation can be 
checked by solving a convex semidefinite optimization problem 
that only depends on the system’s static gain. 

• We illustrate how our method can be applied to study the robust 
stability of the well known Foschini-Miljanic G1 algorithm for 
power control in wireless networks in the presence of uncertain 
interference gains. 

The structure of the paper is the following: in Section|n]we introduce 
the uncertainty setting and the structure singular value, p. In Sec- 
tion|III]we introduce positive and positively dominated systems and in 
Section ITVl we present the results on the structured singular value for 
nonnegative matrices. In Section |V] we derive the convex conditions 
for robust stability of positive systems and finally, in Secion lVll we 
show the applicability of our method to the analysis of robustness of 
the Foschini-Miljanic algorithm. A preliminary version of this work 
was presented in GD, where only positive systems interconnected 
with full block complex uncertainties were taken into consideration. 


Notation 


The set of natural numbers is denoted by N, the set of reals by R 
and K+ (R-I-+) denotes the set of nonnegative (positive) reals. The 
set of complex numbers is denoted by C and Z denotes the integers. 
With ’L[a,b] we denote the set [a, 6] PI Z. The set of Metzler matrices 
(matrices with nonnegative off diagonal elements) is denoted by M. 
The set of nonnegative (positive) diagonal matrices is denoted by D+ 
(]D)++). Given a matrix A G denotes its transpose and 

A* its conjugate transpose. We use d'(A) to indicate the maximum 
singular value of A, tr(A) to denote its trace and p{A) its spectral 
radius. We write A > 0 (A > 0) if A has nonnegative (positive) 
entries and A 0 (A 0) to denote that A is symmetric and 
positive semidefinite (definite) and we denote with |A| the matrix 
of the modulus of the elements of A. For x G C"*, ||a;|| = y/x*x. 
L 2 denotes the space of square Lebesgue integrable signals and Hoo 
denotes the space of causal linear time invariant operators that map 
L^ to itself. For p G the norm is defined as: 



For M G Hoo the operator norm is defined as: 


Halloo 


sup 

pGL^\{0} 


\\Mp\U 

IIPII2 


This norm is sometimes referred to as the Hoo gain of M. We 
denote by M{juj) the transfer function of M in the Laplace domain 
evaluated on the imaginary axis. 



II. Robustness analysis and the framework 

A. Mixed real structured uncertainty 

In this section we introduce the model for uncertain systems that 
will be analyzed throughout the paper. Given /, s G Z[o,oo)i we 
consider two sets of indices I* and Xc such that I* Ulc = 
and XunXc = 0. The uncertainty structure we consider is a set of 
block diagonal complex matrices where the diagonal is composed of 
/ full block uncertainties and s repeated scalar uncertainty blocks. 
The uncertainties corresponding to the index set Ik are real while 
those corresponding to the index set Zc are complex. More formally 
the uncertainty lives in a linear subspace of defined by 

Arc ~ { diag (Ai,..., A/, , ■ ■ ■, ■ 

A. ^ Vi GZRnZ[ij], 

Ai e GZcnZ[ij], (1) 

Si G KVi G Zr n 
Si G CVi G Zc n 

The unit ball in Arc is defined by 

:= {A G Arc : a(A) < 1} . (2) 

Given the preliminary definitions above we consider a linear time 
invariant system M G Hoo, and the following feedback interconnec¬ 
tion: 


p = Mq 
q = Ap, 


(3) 


depicted in Figure [T] where A G Hcxj is an unknown stable linear 
time invariant system such that 


A(jaj) G Sarci ^ (4) 


Note that the blocks in Arc corresponding to Zc represent dynamical 
uncertainties while those corresponding to Zr represent unknown real 
parameters. 

Remark 1: Whenever we consider uncertain interconnection in 
we always assume the system M to be square (the number of outputs 
is the same as the number of inputs). This assumption simplifies the 
notation and is not restrictive as one can pad the transfer matrix 
with zeros to make it square without affecting the stability of the 
interconnection. □ 

We are interested in studying the robust stability of system © with 
respect to all possible stable linear time invariant operators A that 
satisfy the norm bound (0. This is a well studied problem which has 
been proven to be NP-hard in general (H Theorem 2.7]. In this paper 
we consider such problem in tbe special case when M is a positively 
dominated system as defined in Section nni Under this assumption, 
in section |V] we will show necessary and sufficient conditions for the 
robust stability of 0 in the form of Linear Matrix Inequalities. 
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Fig. 1: The interconnection of M and A. 


B. The Structured Singular Value 

Before considering the case where M is restricted to be a positively 
dominated system, we present a brief overview of the classical 
conditions that guarantee the robust stability of 0 in the general 
case as they are useful to develop the main results of the paper. To 
do so we introduce the structured singular value, /i, a widely used tool 
for robustness analysis 03, robust control synthesis 03 and model 
validation da. The framework” provides necessary and sufficient 
conditions for the interconnection in 0 to be robustly stable for all 
stable LTI A satisfying 0. Given a matrix M G and a set 

A C describing the structure of the uncertainty, the structured 

singular value is defined as 03 


^ ■ inf{||A|| : Ag A,det(7-MA) = 0}’ 

if there is no A G A such that det(7—MA) = 0, then p{M, A) = 0 
by convention. The following result due to Doyle [13, GD relates 
the robust stability of the interconnection in 0 to infinitely many p 
tests. 

Proposition 1: Let M G BIoo be a stable LTI operator. Then the 
interconnection in 0 is stable for all A G Sarc if and only if 

sup/i(A7(yu;), Arc) < 1. (5) 

i,;6K 

□ 

Even a single evaluation of /r is a hard problem, according to 
Proposition [T] we need to evaluate p infinitely many times in order to 
obtain a certificate for robust stability. Fortunately upper and lower 
bounds to p are known and they allow to obtain tractable sufficient 
or necessary conditions for robust stability. 


C. Standard bounds of p 

In this subsection we review some important results on the struc¬ 
tured singular value and its bounds. For a more exhaustive review on 
this topic we refer the reader to [T3 or Chapters 7-8]. Given 
an uncertainty structure A, the following upper and lower bounds 
apply; 

p{M) < p{M, A) < (6) 

where © is set of positive definite matrices that commute with all 
A G A defined as 

© := {0 ^ 0 I ©A = Ae, VA G A}. 

The upper and lower bounds in 0 are not tight in general, however 
in Section ITVl we show that if M > 0, then the upper bound is tight 
for the most common block diagonal structures and, in particular, it 
is tight for the structure Arc introduced in 0. 

III. Positive and Positively dominated Systems 
A. Preliminaries on positive and positively dominated systems 

As the goal of the paper is to characterize robust stability for 
positive and positively dominated systems, we now provide the basic 
definitions and the necessary background results concerning these 
classes of systems. 

Definition 1: A linear time invariant operator M is said to be 
externally positive if its impulse response is nonnegative. □ 

It can be shown that a nonnegative impulse response is a necessary 
and sufficient condition to have nonnegative output for all nonnegative 
inputs. External positivity is an input-output property and it is 
independent of the realization. 

We now present some results concerning externally positive sys¬ 
tems. 













Proposition 2 ( hll'i ): Let M{juj) be the transfer function of an 
externally positive system M £ Hoo- Then 

||M||oo = sup(j(M(ja;)) = (j(M(0)), and M(0) > 0. (7) 

□ 

In other words, for externally positive systems, the Hoo gain is 
attained at zero frequency and the transfer matrix evaluated at zero 
frequency is nonnegative. Based on these properties, a larger class 
of systems called positively dominated can be defined. Such systems 
were introduced in El and are defined as follows. 

Definition 2: Let be the element of the transfer 

matrix M{jui) of a system M G Hoo, M is positively dominated if 

sup\mik{ju})\ = mik{0), V(i, fc) £ . ( 8 ) 

□ 

The class of positively dominated systems is strictly larger than that of 
positive systems. Clearly 0 holds for positively dominated systems 
and every positive system is also positively dominated while the 
opposite is not trueQ Note that the definition of positively dominated 
only applies to stable LTI systems, i.e., those mapping L 2 to itself. 
As already noted in (T), at least for rational transfer functions, testing 
whether a system M is positively dominated is a convex problem and 
can be solved efficiently with semidefinite programming. It is also 
interesting to notice that delays, as they do not affect the system’s 
gain, do not change positive domination, that is a system with delays 
is positively dominated if and only if the corresponding delay-free 
system is. 

Remark 2: As already observed in O, the set of positively dom¬ 
inated systems is a convex cone: if G'i(s) and 6*2 (s) are positively 
dominated, then riGi(s) -f T 2 G 2 (s) is positively dominated for all 
'ri,'r2 £ K+. 

B. Exact relaxations for positive quadratic programming 

We now mention a useful theorem which allows us to reformulate 
a class of non-convex quadratic feasibility problems as convex 
semidefinte programs (SDPs) and will be exploited in the next 
section. 

Proposition 3: 1221 Theorem 3.1] Consider Mo, ... ,Mm G M’’ 
and the optimization problems 


1 

f Pt 

= max 

MqX 


Pi ■■ 

ccGR^ 



1 

[ 

s.t. 

MiX > bi 

V i G ^[i,jv] 

1 

{ P2 

= max 

tr(MoA) 


P 2 : 





[ 

s.t. 

ti{MiX) > bi 

V i G Z[i.jv] 


Then the following statements hold true 

i) Pi =P 2 - 

ii) If P 2 has a solution, then it also has a rank 1 solution. 

iii) Let X* be a solution of P 2 , then x = [\/^Tli • ■ • > s/^NN ] 

is a solution of Pi. □ 

IV. Structured singular value for 

NONNEGATIVE MATRICES 

In this section we consider the structured singular value for 
nonnegative matrices. In particular we show that if M > 0, then the 
inequality 0 always holds with equality for the upper bound. The 
structured singular value p can therefore be computed by solving an 

*The system ^ with -M < C < 1 is positively dominated 

+ \/2 f J 

but not externally positive since the impulse response changes sign infinitely 
many times. 


LMI problem. Before proving such equality we need the following 
result: 

Theorem 4: Let M £ and Arc defined in 0 and 

s := If /i(M, Arc) > 1 then there exist a A = 

diag(Ai,..., A/, 5/s) £ AucnR)))^"*, with ||A|| < 1 such that 
det(/ — MA) — 0, furthermore there exist a nonnegative vector 
q £ Rip such that q = MAq. □ 

Proof: By Assumption p{M, Arc) > 1, therefore there exists 
A £ Arc such that det(/— MA) = 0, and ||A|| < 1. We conclude 
that p{MA) > 1. We now consider the following matrix: 

A :=diag(|Ai|,...,|A/|,57s), 

Where S := maxj |5j|. Notice that A £ Arc DR™^™. From 1231 
Theorem 8.1.18] we know that: 

1 < p(MA) < p(|MA|) < p(M|A|) < p(MA). (9) 
Now we need to study the norm of A. Since A £ A we know that: 

A = diag(Ai ,..., Af, , ■ ■ ■, Sf+slmf^^), 

and ||A|| = max{maxfc ||Afc|| ,maxj |5j|} < 1. It is well known 
that every block At can be taken as a dyad (matrix of rank 1), see 
for example |19| . We therefore write Ak = ^kCk ’ clearly ||Afc|| = 
ll6IIIICfe|| = III6III IllCfclll = ll|Afe|||. We conclude that: ||A|| = 
IIAII < 1. From 0 we notice that A := p(MA) > 1. We now define 
the matrix A := A/A, clearly || A|| < 1 and p(MA) = 1. Since MA 
is a nonnegative matrix, we know from Perron-Frobenius Theorem 
that there exist a nonnegative eigenvector q such that MAq = q, 
which implies that {I — MA)q = 0. Therefore {1 — MA) is singular 
and and the proof is complete. ■ 

Theorem |4] has two important implications, which we summarize 
in Corollary [3 

• It allows us to replace possibly complex valued scalar blocks 
with full blocks in the perturbation A (with some blocks of 
dimension 1). 

• It allows to consider real nonnegative perturbations without loss 
of generality. 

Corollary 5: For M > 0 robust stability with respect to the 
uncertain set: Bakc defined in 0, is equivalent to robust stability 
with respect to the set: 

Z3a = {diag(Ai,...,A/,5/+i,...,5/+s), 

Afe £Rp‘'=^"‘'=VA:£Z[ij], (10) 

Sk G R+VA: £ , ||A|| < 1}. 

where s := X)fci/+iOT in other words, p(M, Arc) = 
p(M,A). □ 

Proof: By Theorem [T] robust stability with respect to Ba and 
robust stability with respect to Bakc ^re both equivalent to robust 
stability with respect to 

^A { diag(Ai,..., Af,5Is), 

A^ £Rm,xmfcVA:£Z[ij],5£R+, ||A|| < 1}. 

therefore /i(M, Arc) < 1 <=4- p(M, A) < 1. Since p{aM,-) = 
\a\p{M,-) m, it is clear that, unless p(M, Arc) = p(M, A), 
we could find a value a > 0 such that p{aM, Arc) = 1 and 
p{aM, A) < 1 reaching a contradiction. ■ 

A. Block diagonal uncertainty structure 

In view of Corollary |5] when dealing with nonnegative matrices, 
we can restrict ourselves to study robust stability with respect to the 
block diagonal real-nonnegative uncertainty structure 

A := {diag(Ai,..., An) : Ak G Vfc £ Zfi^jyj}, 




which contains only full-blocks (possibly of dimension mk=l for 
some fc), as all cases including scalar-times-identity or complex 
blocks can be re-casted in this form w.l.o.g. 

We define the set © of all positive definite matrices that commute 
with A given by 

© — {diag(6li7,„i,... ,6'jv7m„) : £ K++VA: G 

We now need the following lemma: 

Lemma 6: Let p,q & IR+. Then there exist a matrix A £ 
with IIAII < 1, such that q = Ap if and only if Ikll < IIpII- o 
Proof: We assume that ||q|| < ||p||, then 



satisfies the requirements. The other direction comes directly from 
the definition of matrix norm. ■ 

We can use Lemma to characterize necessary and sufficient 
conditions for p{M, A) < 1 in terms of inequalities in the norms 
of the components of q and Mq corresponding to the blocks in the 
structure of A. We can write these components as Ekq and EkMq 
respectively, where 

Ek ■= [Omi • ■ • Irrif, ]. 

Using this, we can define the quadratic functions 

Ml) ~ UMq)kf - \\qkf = q^ (M^EJEt:M - EJEt:)q, 
and the set 

:= {{(j>i{q),.. . ,</>jv(q)) : q £ R+, ||(j|| = 1} • 

With the following Proposition we can give a geometrical character¬ 
ization for the robust stability test. 

Proposition 7: For any M £ The following statements 

are equivalent. 

(a) p{M, A) < 1, 

(b) The sets Vr^ does not intersect the positive orthant. 

□ 

Proof: We will show that -■(a) The set Vr.,. 

intersects the positive orthant if and only if there exists q £ R+, 
with ||g|| = 1 such that: fkiq) > 0, Vfc £ jv]- By Lemma|^ this 
happens if and only if, Vfc £ 'Z[i,n] there exist At, with || Afc|| < 1, 
such that 

AkEkMq = Ekq, V fc £ 

Putting all these blocks together this is equivalent to the existence of 
q, with ||q|| = 1 and A £ A with ||A|| < 1 such that: AMq — q, 
which is equivalent to the fact that {I— AM) is singular, which is in 
turn equivalent to the singularity of (7 —A7A), which is the negation 
of (a). ■ 

Lemma 6 and Proposition 0 are the real nonnegative counterpart of 
known results that holds for complex valued M and A 1201 Lemma 
8.24, Proposition 8.25], 1241 . l25l . The fundamental difference is 
that, if M £ we can w.l.o.g. restrict the search for q from 

C"* to R+. This enables us to exploit powerful tools from nonlinear 
optimization to find tractable conditions for robust stability leading 
to the main result of this section. 

Theorem 8: Let M £ R(|]^"‘. Then the following are equivalent: 

a) p{M, A) < 1 

b) There exists 0 £ © such that ||02M0“5 || < 1. 

□ 

Proof: We define the matrices := E^EkM—E^Ek. In 
Proposition|7]we established that (a) is equivalent to the fact that Ve_|_ 
and the positive outhant are disjoint. We can re-write this condition 


in terms of the infeasibility of the following non-convex quadratic 
program 

q^ Miq>Q, Vi£Z[i_jv] 

q^q = l ( 11 ) 

q £ R+. 

We notice that the since TJj Ek is a partition of the identity and thus 
diagonal, {Mk} are Metzler matrices. By Proposition 0 we know 
that the standard SDP relaxation for non-convex quadratic programs 
is exact. We can therefore reformulate the program in GD as 

tr(MiQ) > 0, Vi£Z[i_]v] 

tr(g) = 1 (12) 

Q 0, 

and if Q = \qij] solves (T^, then q = [y'qTT,..., ^qNN] solves CB- 
Infeasibility of is thus equivalent to infeasibility of We can 
now use [Bl Proposition B.2], which exploits a version of Farkas 
Lemma for LMIs, to prove that J12l l is infeasible if and only if there 
exist multipliers > 0 such that 

N 

Y^OkMk^Q. (13) 

i=k 

Since the matrices Mt are nonnegative except possibly for a partition 
of the diagonal corresponding to —EjEk, the matrices EJEk 
together form a disjoint partition of the identity and a necessary 
condition for dkMk ^ 0 is that the whole diagonal is negative, 

the multipliers {6k} must all be strictly positive for ( I13t to hold. We 
then, w.l.o.g., restrict the search to 9k > 0. The condition in l ll3t can 
thus be rewritten as follows. There exist dk > 0 such that 

N 

Y,ekM^EjEkM-ekEjEk-<0. (14) 

i = k 

Since 

N 

OkEjEk = diag(0i7,..., 0jv7) £ ©, 

i = k 

the condition in J14l l can be equivalently reformulated as 3 0 £ © 
such that 

M^0M-0^O. (15) 

From simple algebraic manipulations one can show that the condition 
in J15t is equivalent to the existence 0 £ © such that 

||05M0-i|| < 1. 

which concludes the proof. ■ 

One interpretation of the multipliers {6i} is that of the coefficients 
of a hyperplane that separates the sets Vr_,_ and the positive orthant. 
From the proof we can deduce that, if M > 0 despite the fact that 
Vr^ is not necessarily convex, if p,{M, A) < 1, then there exist an 
hyperplane that separates Vr^ from R^|^. 

The following corollary is an immediate consequence of Theo¬ 
rem!] 

Corollary 9: For all M £ R(|)^"*: 

p(M, A) = Jnf^||0^A70-^||. (16) 

Since p{aM, A) = \a\p{M, A) II19I . if | |161 did not hold one could 
scale M in a way that contradicts Theorem!) 

Notice that the inequality GB is an LMI and therefore it can 
be solved in polynomial time using an interior point algorithm. An 
early result that shows that the convex upper bound to p is tight for 
nonnegative matrices appears in 1261 for the particular case when the 
matrix A is restricted to be diagonal. 




V. Robust stability of positively 

DOMINATED SYSTEMS 

All results so far were dealing with matrices only, now we want 
to obtain conditions for robust stability of positively dominated and 
positive systems. Let M € Hoc, according to Proposition [T] in order 
to test the robust stability of the interconnection in 13 with respect of 
all A € Bakc’ rieed to check that sup^^gj, Amc) < 1- 

We propose a frequency domain characterization of robust stability 
for positively dominated systems and we show that we can obtain a 
convex condition involving only the transfer matrix at zero frequency. 
First we need to define a new uncertainty set related to Arc- Given 
Arc defined in m and s consider the sets 

A = { diag(Ai, - - -, A/, Sy+i,6f+s) ■ 

Ak Gvfe G 
5k G K+Vfc € 

and the set of positive definite matrices commuting with A 

G) = { diag(0i/7ni 5 • • ■) ^ ■ ■ ■ 5 ^/+s) ■ 

Ok G 

We already showed, with Corollary that for the case of positive 
matrices, robust stability with respect of the unit ball of the set Arc 
is equivalent to robust stability with respect of the unit ball of the 
set A. In the following we show that this still holds for the case of 
positively dominated systems. 

A. Frequency domain convex characterization 

The following theorem is the main result of the paper and provides 
a convex characterization for robust stability of positively dominated 
systems in the frequency domain. 

Theorem 10: For a positively dominated system M £ Hoo then 
the following holds: 


sup^(M(jtj), Arc) = sup/r(M(jtj), A) = inf ||02M(O)0 ^ ||. 

I.J6M i-jSR 

Proof: We start by considering ^(M(0 ),Arc). Since M is 
positively dominated, we have that M(0) £ Therefore, 

we can apply Corollary |5] to consider full block real nonnegative 
uncertainties only and Corollary[9]to show equivalence with the upper 
bound and we get 


/i(M(0),Arc) = /r(M(0),A) = ||02 M(O)0-5 | 


(17) 


This implies that there exist a sequence {0*, £ ©}fegN such that 
lim ||0 |m(O) 0“5|1 =/r(M(0),ARc) =ai(M(0),A). (18) 

k—^co 

Using the upper bound given by at each frequency, we have the 
following inequalitieas for all 0 ^ in the sequence 


sup/i(M(yu;), Arc) < sup ||0 = M(ja;) 0 j, =||, and 


cjGR 




sup jut), A) < sup||0^M(jcu)0j,' 


(19) 


Since 0| £ © C D++ and set of positively dominated systems is 

a cone, 0 |M(ya;) 0 j, ^ is itself the transfer function of a positively 
dominated system. From Proposition [3 we have 


sup||0|M(ja.)0,^|l = ||0|M(O)0“' 


( 20 ) 


We now let fc —^ oo and, in view of Cl, ODi and (|23 we have 

sup/i(M(ja;), Arc) < ^(M( 0 ), A) = ^(M(0), Arc), and 

. ( 21 ) 

sup^i(M(ja;), A) < A) = /i(M(0), Arc). 

The inequalities in (EB always hold with equality as 
m(('^(0), Arc) < sup„gR^i(M(jcu), Arc) and ^i(M( 0 ),A) < 
supj^gjf/i(M(jtj), A) respectively. In view of l ll7t the proof is 
complete. ■ 

In view of Proposition [T] Theorem [TO] gives us a convex character¬ 
ization of robust stability for uncertain positively dominated linear 
systems with respect to all perturbations in -Barc- Note that the 
positive-domination assumption is required only for M and not for 
the whole interconnection in l| 3 - 

VI. Application to the robust stability of the 
Foschini-Miljanic algorithm 

We apply our result to characterize the robust stability of the widely 
studied distributed power control algorithm for wireless communica¬ 
tion presented in m by Foschini and Miljanic. 


A. Review of the Foschini-Miljanic algorithm 

We consider N communication channels, each corresponding to a 
reciever-transmitter pair. We denote by {hi} £ (0,1] the transmission 
gain for channel i and by {gji}iy:j £ [ 0 , 1 ] the interference gain from 
channel j to channel i. We denote by Oi the variance of the thermal 
noise at receiver i and by pi the power level chosen by transmitter i. 
The quality of the transmission at channel i is measured by the Signal- 
to-Noise-and-Interference-Ratio (SINR), denoted as Fi and defined 
as 

p _ hiPi 

9jiPj + n' 

Assuming that each receiver can measure its SINR F i and communi¬ 
cate it to the transmitter, the aim of the Foschini-Miljanic algorithm 
is to chose the power levels {pi}{Li such that Fi > ji for all 
i £ Z[i jv], where 7 ^ £ R+ are pre-defined constants that guarantee 
a certain quality-of-service for each communication link. 

The algorithm reads as follows: 


Pi = ki 



9jiPj + I'i 




( 22 ) 


where ki £ R++. The algorithm is decentralized in the sense that 
each transmitter has access only to the interference measurement 
of the receiver: F := 'Y^jy.iPjiPj + Fi. If we define the 
vectors p = (pi,...,pjv)^, v = (ui,..., Ujv)^, and the matrices 
K = diag(fei,..., fciv), (& = diag(^,.. ., ^) and G such that 


Gi ■ = I ° ‘f * = 7 

\ gji otherwise ’ 

the algorithm in d 22 b can be rewritten in matrix form as 


p = K{-I + T'G)p -b K'Fu, (23) 


Since K{—I + T>G) £ M and K'^ > 0, (|23 describes a 
positive system. In 03 it was proven that, if there exists a vector 
p = (pi,...,P jv)^ such that 


hiPi 

T,jjii 9jiPj + n 


> 7i 


Vi £ Z[i^jv], 


(24) 


Then the algorithm in l l22b will converge to such a solution. If a 
solution p satisfying El does not exist, then E3 will be unstable. 




B. Robust stability of the Foschini-Miljanic algorithm 

We now show that the analysis tools developed in Section lYl can 
be used to establish whether the Foschini-Miljanic algorithm (FM 
algorithm) is robustly stable when the interference matrix G is not 
known exactly but can be described in a parametric form 


G = Go + SAF, (25) 

where A has an arbitrary block diagonal structure and (f(A) < 1. 
We can consider the system 


r K{-I + <^GQ)p + K<^Eq 

[ z = Fp 


and the interconnection 


p = Mq 
q = Ap. 


(27) 


Since M is positive and thus positively dominated, if we consider 
M(0) = F{—I + ^'Go)“^5'i? to be its transfer function transform 
evaluated at zero frequency, we know from Theorem 1 101 that we have 
robust stability with respect to all uncertainties satisfying l l25t if and 
only if 30 € © such that 

M(0)^eM(0) -0^0. ( 28 ) 

Note that, with standard p analysis, the LMI test in J28t would be 
neither necessary nor sufficient to establish robust stability of the 
FM algorithm with respect to the uncertainty in the matrix G. Not 
necessary because the p upper bound is in general not tight and not 
sufficient because the test is performed only at zero frequency. In 
contrast with standard p analysis, by exploiting the positivity of M, 
and applying Theorem 1 101 we know that the LMI test in ( I28t is both 
necessary and sufficient to test robust stability of the FM algorithm 
with respect to the uncertain interference matrix G. 

In the literature it was shown that the delayed FM algorithm 
is stable for any bounded delay if and only if the delay-free system is 
stable. A similar result can be easily established for robust stability 
with respect to the interference matrix G using the analysis tools 
developed in this paper. Adding delay elements of the form to 
the transfer matrix preserves the positive domination of the system 
and, at the same time, does not change the zero frequency matrix 
M(0). The stability test remains therefore unchanged in the presence 
of delays implying that the delayed FM algorithm is robuslty stable 
with respect to G if and only if the corresponding non-delayed 
algorithm is robustly stable. 


VII. Conclusion AND OUTLOOK 

In this paper we provide tractable necessary and sufficient condi¬ 
tions for robust stability of uncertain LTI systems modeled as the 
feedback interconnection of a positively dominated system and a 
block structured, norm bounded stable LTI system. We give LMI 
conditions that only depend on the static gain matrix. Our results 
falls in line with several recent results on positive systems showing 
that problems that are hard for general systems become easy when 
restricting to this simpler class of systems. The results we propose are 
expressed in terms of LMIs; this is in contrast with most of the recent 
literature on positive systems, where results are obtained prevalently 
in the form of Liner Programs (LP), see for example m, (a, m- 
LPs allow greater scalability as they can be solved for larger systems 
and to higher accuracy. An LP reformulation of our robust stability 
test is an interesting problem to investigate as it would allow scalable 
robust stability verification for very large scale systems and it is left 
as future research. 
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